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A metric-afRne approach is employed to study higher-dimensional modified gravity theories in¬ 
volving different powers and contractions of the Ricci tensor. It is shown that the field equations are 
always second-order, as opposed to the standard metric approach, where this is only achieved for 
Lagrangians of the Lovelock type. We point out that this property might have relevant implications 
for the AdS/CFT correspondence in black hole scenarios. We illustrate these aspects by considering 
the case of Born-Infeld gravity in d dimensions, where we work out exact solutions for electrovacuum 
configurations. Our results put forward that black hole singularities in arbitrary dimensions can be 
cured in a purely classical geometric scenario governed by second-order field equations. 
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I. INTRODUCTION 

A major reason to consider the physics of modified 
gravities in extra dimensions is provided by the Anti-de- 
Sitter/Conformal Field Theory (AdS/CFT) correspon¬ 
dence. This is a duality relating the physics (more pre¬ 
cisely, the thermodynamics) of black holes defined in AdS 
space (n-dimensional) and the CFT lying on its ((n— 1)- 
dimensional) boundary [l|. It constitutes an example 
of the more general concept of gauge/gravity dualities 
and holography, by which the physics of the interior of 
a system can be described through the physics of its 
boundary. Even though the AdS/CFT correspondence 
was originally developed within string theory, it has far- 
reaching consequences and applications in other areas, 
such as in the description of strongly-interacting and con¬ 
densed matter systems 0 ■ One can thus forget about the 
original motivation for the correspondence, focus on the 
analysis of properties of black holes in AdS space-times, 
and try to work out the corresponding aspects of the CFT 
using the dictionaries established in the literature Q- 

On the CFT side, one of the most interesting cases 
is the four-dimensional one, implying that one is led to 
consider gravity theories in five dimensions. An impor¬ 
tant lesson of the existing AdS/CFT dictionaries, is that 
new couplings and interactions in the gravitational sec¬ 
tor added to the Einstein-Hilbert (EH) Lagrangian of 
General Relativity (GR) have a direct correspondence on 
the CFT side. By adding new elements in the gravita¬ 
tional sector one could thus explore new couplings in the 
CFTs. However, the addition of higher-order curvature 
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invariants in the gravitational sector usually gives rise to 
higher-order derivative field equations. This largely pre¬ 
vents finding analytical solutions and usually results in 
the appearance of ghosts and instabilities, which casts 
severe doubts on the physical consistency of the corre¬ 
sponding theory. 

A well known way out of this problem is to consider 
a particular combination of quadratic corrections to the 
EH Lagrangian that removes the undesired higher-order 
derivative terms, known as Gauss-Bonnet theory. In this 
theory, in addition, exact solutions can be obtained [^. 
This theory contains a single independent parameter, 
which limits the range of dual CFTs that can be studied. 
A natural next step is to add cubic curvature terms to 
the Gauss-Bonnet theory in appropriate combinations so 
as to keep second-order field equations. This procedure 
can be generalized to arbitrary orders by means of the 
so-called Lovelock theories A drawback of using the 
Lovelock approach to add cubic couplings is that in five 
and six space-time dimensions the corresponding cubic 
terms become topological invariants. This means that 
they add no new contribution to the field equations (the 
same applies for the Gauss-Bonnet combination in four 
space-time dimensions). One is thus faced with a seri¬ 
ous limitation to include new curvature couplings in the 
gravitational sector which can contribute in a nontriv¬ 
ial manner to the field equations, which would have an 
impact on the dual CFT. 

To circumvent this argument, a new procedure has 
been developed based on the addition of a set of new 
curvature-cubed terms to the action in five dimensions 
such that, though the field equations are, in general, of 
third-order, the linearized equations describing gravitons 
propagating in the AdS background are second-order . 
They are called “quasi-topological gravities”, as the extra 
terms are not topological invariants any longer. How¬ 
ever, besides the lack of naturalness in this procedure, 
as a matter of fact one is forced to solve highly compli- 
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cated equations where solutions are only obtained under 
restrictive conditions. Extra freedom in these theories is 
thus limited and a aesthetically more appealing approach 
would be desirable. 

In this work we consider a different approach to this 
problem and face it from the perspective of metric- 
affine (or Palatini) geometry. In the Lovelock and quasi- 
topological approaches, one implicitly assumes the hy¬ 
pothesis that the connection is metric-compatible, i.e., 
that it is defined a priori in terms of the Christoffel sym¬ 
bols of the metric. This assumption, which is more a 
historical convention than an experimentally supported 
fac10, is at the root of the limitations of the previous 
approaches. We will see that relaxing the compatibil¬ 
ity condition between metric and connection and allow¬ 
ing the connection to be determined by the field equa¬ 
tions generically produces ghost-free, second-order metric 
field equations. This strategy, therefore, provides a fair 
amount of freedom in the gravitational couplings without 
the restrictions of Lovelock theories and quasi-topological 
models. 

It is worth stressing at this point that metric and con¬ 
nection carry very different physical meanings. While 
the former deals with properties locally defined in space, 
such as measurements of lengths, areas, volumes, etc, the 
latter is related to properties remaining invariant under 
affine transformations, such as parallelism (see for a 
pedagogical discussion). Forcing these two structures to 
be related to each other a priori is an unnecessary con¬ 
straint that limits the potential of these theories. In the 
metric-affine viewpoint, the field equations are obtained 
by independent variation with respect to metric and con¬ 
nection. As we will see, in many cases of interest, the 
connection turns out to be a constrained (non-dynamical) 
object which, in general, differs from the Levi-Civita con¬ 
nection (the exception to this rule corresponds precisely 
to Lovelock theories ©[HI). As GR is a particular 
case of Lovelock gravities, one finds that the Palatini ap¬ 
proach gives in this case the same results as the standard 
metric approach. However, in general, the mathematical 
and physical properties of modified gravities in the met¬ 
ric and Palatini approach largely differ from each other 
(see [l^l for a review on the topic). 

As already mentioned, a key aspect of the Palatini for¬ 
mulation is that, for a large family of functional forms of 
the Lagrangian density, the field equations turn out to 
be second-order. This is so due to the active role played 
by the matter in the construction of the connection. In 
particular, in absence of matter fields, the connection 
becomes metric compatible and the field equations boil 
down to those of GR with a cosmological constant term. 


^ Indeed, analysis of the physics of crystalline structures with 
defects on their microstructure, supports the view that non- 
Riemannian geometries are favoured in Nature Q, which might 
have important consequences for the understanding on the mi¬ 
croscopic structure of space-time @]. 


which avoids the presence of extra propagating degrees of 
freedom. When matter is present, the connection equa¬ 
tion can be seen as defining the Levi-Givita connection 
of a new metric , which is algebraically related to the 
metric originally appearing in the action. This al¬ 
gebraic relation depends on the matter sources and on 
the particular gravity Lagrangian chosen. Additionally, 
one finds that the field equations for hp^ can be cast in 
Einstein-like form, which greatly simplifies their resolu¬ 
tion. 

Recently, we have successfully applied this formalism 
to the case of higher-dimensional f{R) Palatini theories 
and worked out exact solutions in the five-dimensional 
case [l^. Now we extend the scope of the methods de¬ 
veloped there by considering the case of theories includ¬ 
ing also powers of the Ricci-squared scalar, Q = Rp^R^'', 
and other curvature contractions. In particular, we con¬ 
sider a gravitational Lagrangian of the Born-Infeld type, 
which formally contains up to four powers of the Ricci 
tensor. We derive the field equations for these theories 
and write them in terms of the hpv metric to highlight 
their second-order character and how a cosmological con¬ 
stant naturally emerges. Focusing on Born-Infeld gravity, 
we work out exact solutions corresponding to an elec¬ 
trostatic spherically symmetric electromagnetic field. A 
discussion on the internal structure of its solutions and 
its implications is provided. In particular, we show that 
these solutions are self-gravitating structures with non¬ 
trivial (wormhole) topology and, at the same time, are 
geodesically complete and non-singular. This is so de¬ 
spite the presence of curvature divergences, which de¬ 
mands for a reconsideration of the role traditionally at¬ 
tributed to curvature invariants to characterize space- 
time singularities, since their presence in our scenario 
pose no obstacle for the extendibility of paths through 
the wormhole. 

Regarding this last point, we note that one’s intuition 
tends to establish a correlation between space-time singu¬ 
larities and curvature divergences. In fact, this question 
has been discussed in the literature from different per¬ 
spectives, including the viewpoint from the philosophy of 
science (see 0 and references therein). However, from a 
mathematical perspective, the characterization of singu¬ 
lar space-times is more subtle [l5l - 0 . In this sense, even 
Minkowski space-time with a point (or small domain) re¬ 
moved can be regarded as singular despite the complete 
absence of curvature divergences. The key point is that 
the removal of a point (domain) has a deep impact on the 
existence of physical observers. If observers are regarded 
as point-like entities that follow geodesics, the removal of 
a point (domain) implies that observers reaching to that 
point (domain) may simply cease to exist there (or come 
into existence from nowhere). Since that point (domain) 
does not belong to the space, physical observers are sim¬ 
ply not defined there. There is no Physics at that point 
(or domain). Whether an extended object can experience 
arbitrarily high accelerations or deformations due to cur¬ 
vature divergences is of secondary importance as long as 
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it can exist. The key point, therefore, is the existence of 
physical observers, which are characterized by complete 
geodesics. 

In a series of previous papers, some of us addressed the 
issues of geodesic completeness in black hole space-times 
with wormhole structure the impact of curvature 
divergences on extended observers (defined by geodesic 
congruences) [l9j| . and also the characterization of cur¬ 
vature divergences by means of the scattering of (scalar) 
waves off the effective potentials that they generate . 
It was found that these space-times are geodesically com¬ 
plete, that the volume of extended objects remains finite 
as the divergent region is crossed, and that wave prop¬ 
agation is uniquely determined across the curvature di¬ 
vergence, being possible the computation of the reflec¬ 
tion and transmission coefficients that characterize the 
barrier associated to the curvature divergence. In this 
work we will extend some of those results derived in four 
space-time dimensions to the higher dimensional case. 


II. PALATINI THEORIES WITH RICCI 
SCALAR AND RICCI-SQUARED INVARIANTS 

In order to illustrate the basic strategy to deal with 
Palatini theories, consider d-dimensional theories defined 
by the following action 

5=^/ (1) 

where is a constant related to d-dimensional New¬ 
ton’s constant in appropriate system of units, g is the 
determinant of the space-time metric the curvature 
invariants R = gfj.^R^'' and Q = R^^R^^'' are constructed 
with the Ricci tensor, R^^ = R^^pv^ which follows from 
the Riemann tensor 




p.0 
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In these expressions the connection T = is a priori 
independent of the metric g^^. For simplicity, we as¬ 
sume that the connection is not coupled to the matter 
action, Sm, which depends only on the metric and on 
the matter fields, denoted collectively as tpm- In general, 
the antisymmetric parts of both the connection (called 
torsion), and the Ricci tensor, R[p,j], can be non¬ 

vanishing. In this work, however, both contributions are 
neglected for simplicity (see [2lj for a discussion of the 
general case). Under these conditions, variation of the 
action © with respect to metric and connection yields 


<5S'=^ J {fRRpu + ‘^fQ9°'^RpaRpu- ]^griuf)5g^^ 


-{fRg^'^ + 2fQR^pg<^>^gP'')5Rp,{T) 


SS„ 


( 3 ) 


where we have used the short-hand notation fx = 
df/dX. An important point now is to realize that the 


Ricci tensor is an object constructed only with the con¬ 
nection and, therefore, does not depend on the metric g^^. 
This point is crucial to understand the different structure 
of the field equations as compared to the (more standard) 
metric approach. In the Palatini formulation, we have to 
write the variation of Rpu in terms of the variation of 
the connection. To do this we make use of the relatiorU 
dRpuiX) = (valid for torsionless con¬ 

nections [13), and consider the piece of the variation in 

as 

= ^ / d'^x^gMp.iyxSTl^ - V.JT^^) (4) 


where we have introduced the object 

= fRg^’^ + 2fQg^^gf^’'R^p (P) (5) 


Integrating by parts we obtain 
1 


6rS = 


d<^x^g 


2k2 

+ + Lsy p{^gM''P) <5r^^6) 


I 


Contracting A and v in the expression above leads to 
(1 — d)V\{y/—gMP‘^) = 0. In summary, the system of 
metric and connection equations for the action © in the 
Palatini formalism with the assumptions above becomes 

IrRpu + 2fQg^^Rp^Rp, - ^gp,f = R^Tp, (7) 

Va(v/=^M'^")=0, (8) 

where Tpi, = —is the energy-momentum tensor 
of the matter. 

To solve these equations we must first show that 
is linear in the connection. At first sight this is not ob¬ 
vious and, in fact, the definition ([5]) suggests that © is 
a nonlinear equation with up to second-order derivatives 
of the connection. The point is that one can use ([7]) to 
show that the object Pp‘' = Rpa{T)g°^'^ is an algebraic 
function oiTp'^, which allows to get rid of the connection 
dependence of in favor of a dependence on the mat¬ 
ter fields and the metric. To see this, we rise one index 
in 0 with gP'^ to obtain 


fRPp'' + 2/QP^“Pa" - = k'T/ . (9) 

By noticing that R = Pp^ and Q = Pp'^P^^, it follows 
that Pp‘' must be an algebraic function of Tp''. Eq.0 
can thus be written as 


VA(\/=^5'^“Sa")=0, (10) 


2 For simplicity we omit the contribution of torsion in this expres- 
sion because we set it to zero at the end. However, one should 
keep those terms in the variation to consistently obtain the field 
equations. See for details. 
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where 


Sa" = IrSc’' + 2/qP„" (11) 

can be regarded as a function of the matter fields, i.e., 
P = P{T), with a hat denoting matrix representation 
of Pa‘^ and Ta‘^. In this form, the connection in m 
appears linearly and can be solved by means of elemen¬ 
tary algebraic manipulations, in much the same way as 
the equation ^= 0 is solved in the Palatini 
version of GR (see chapter 21 of for details). This 
can be seen by just proposing the existence of a rank-two 
tensor such that 


, ( 12 ) 

which turns (1101) into V\{\/—hh^'') = 0. This shows that 
the independent connection can be written as the 
Levi-Civita connection of the auxiliary metric Note 

that this formal solution is valid for arbitrary Lagrangian 
/(i?,Q). 

The explicit relation between and can be ob¬ 
tained by noting that h = 51 S | , where | E | represents 
the determinant of the matrix E. Replacing this result 
back in m one obtains 

h=\t\^±-^g-,h-^ = ^r^ (13) 

|E|^ 

The metric field equations © can be conveniently writ¬ 
ten in terms of by noting that (ED can be written 
as 

P^“E„" = . (14) 

Using the relations dH, one easily verifies that 

P^“E„" = i?^„(r)/i“'^|E|^ = + Given 

that Rfj,a{T) — Rfj,a{h), we finally get 

= . (15) 

where Rfi^{h) = Rf^ah°‘'^ . 

Several comments are in order. First we note that this 
formal representation of the field equations recovers GR 
as a particular case when f{R,Q) = i?, because in that 
case E = / and 5 ^ 1 , = On the other hand, it puts 
forward that the metric satisfies a system of second- 
order field equations with the matter sources on the right 
hand side (recall that R and Q are functions of T'^'^). 
Since the physical metric is algebraically related to 
via the matter-induced deformation , it follows 
that is also governed by second-order equations. To 
explore the properties of vacuum solutions, it is useful to 
complete the square in ([^ and put it (in matrix notation, 
for simplicity) as 






1 + 


2/q" 


-T . 


(16) 


In vacuum, T = 0, it follows that P = a{R, Q)i, where 
the explicit form of a{R, Q) is not relevant for the dis¬ 
cussion. From this relation one finds two constraints, 
R = Tr[P] = 4a{R,Q) and Q = Tr[P^] = 4a(R,Q)^, 
which imply that in vacuum R and Q must be constants. 
We thus find that E^*^ oc Sfj,'', with the proportionality 
factor being a constant, which allows to write m as 
Rfiuig) = ^effgiiv This shows that regardless of the 
/(P, Q) Lagrangian, the vacuum equations of the the¬ 
ory exactly recover the GR equations with an effective 
cosmological constant Ag//, whose form depends on the 
particular /(P, Q) function chosen. 

We have just seen that in the Palatini framework 
/(P, Q) theories provide second-order modified dynamics 
with the same number of propagating degrees of freedom 
in vacuum as GR. The presence of matter fields induces 
a nontrivial deformation between the metrics 5 ^ 1 / and 
which can have remarkable consequences in black 
hole scenarios, such as the replacement of the central 
singularity by a wormhole [l3| . The second-order char¬ 
acter of the field equations is guaranteed for arbitrary 
/(P, Q) Lagrangian, polynomial or not, which contrasts 
with the tight constraints that appear in the usual metric 
formalism and motivate Lovelock theories and the quasi- 
topological models mentioned in the introduction. One 
is thus free to introduce new curvature couplings in the 
gravitational action without worries about the genera¬ 
tion of ghosts or higher-order derivatives. The price to 
pay, however, is the need to solve algebraic equations 
required to express P and Q as functions of the mat¬ 
ter fields and the coupling parameters. Some analyti¬ 
cally tractable models of the /(P) type have already been 
discussed in [T^. The introduction of a Q—dependence 
makes the analysis a bit harder, though some tractable 
ad hoc models can be found. In the next section we con¬ 
sider a better motivated theory which admits a complete 
analytical treatment and possesses all the algebraic prop¬ 
erties found in the /(P, <5) models presented here. 


III. BORN-INFELD GRAVITY IN d 
DIMENSIONS 


An interesting proposal for a high-energy gravitational 
action that recovers GR at low curvatures is that of 
Born-Infeld gravity. The idea is to construct a gravi¬ 
tational analog of the non-linear (classical) extension of 
Maxwell electrodynamics introduced by Born and Infeld 
[^ . where the Maxwell action is replaced by a square- 
root form. In this way, both the divergence of the field 
strength tensor and the self-energy associated to charged 
particles are removed within a classical model of the elec¬ 
tromagnetic field. Deser and Gibbons [ 2 ^ considered 
its gravitational counterpart by essentially replacing the 
field strength tensor by the symmetric Ricci tensor ob- 
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taining 


S = -^ J d'^x + ei?^i.(r)| - \y/^ 

m) 


(17) 


where the same notation as in the previous section ap¬ 
plies. The constant e is a characteristic parameter of 
the theory with dimensions of length squared, and whose 
sign and value depend both on theoretical grounds and 
on its compatibility with experiments. This theory has 
attracted much interest in the last few years, with ap¬ 
plications in astrophysics [ 1 ^, cosmology [l^, and black 
hole scenarios [2^. This theory has also motivated new 
families of high-energy modifications of GR [1^ . 

The meaning of the parameter A in dm) follows from 
an expansion of the gravity sector, Sbi, in series of e —>■ 0 
as 

lim Sbi = ^ J d‘^xy/^[R - 2Aeff] (18) 

■ i ("T + +■■■ 

where A^// = plays the role of the effective con¬ 
stant of the theory. This series expansion allows to see 
that BI gravity recovers the action of GR plus a cosmo¬ 
logical constant in d dimensions at low curvatures, with 
the next-to-leading order corresponding to that of a spe¬ 
cific quadratic f{R,Q) theory. Higher orders in this ex¬ 
pansion involve contractions and products of increasing 
powers of the Ricci tensor. 

The Lagrangian density appearing in m can be con¬ 
veniently rewritten by defining = g^a (<5“y -b eP°‘,j), 
where = g°‘^Rpv{T), and denoting = g^^^qpu = 
i5“iy -|-eP“j,. With this definitions, we have H = g~^q and 
= q~^g (where a hat denotes matrix representation), 
which turns (HZl) into 


Sbi = 


d'^Xy/^ 




-A 


+ Sn 


(19) 


This also suggests that BI gravity is a particular exam¬ 
ple of a more general family of theories, defined by the 

Lagrangian density Lq = /(|^|), with / = the BI 

case (see Ref. [30l|l. 

Following the same procedure as in the f{R, Q) case, 
variation with respect to metric and connection leads to 
the following field equations 

( 20 ) 

yx{^f^qqn = 0 . ( 21 ) 

From the connection equation one readily sees that the 
connection is compatible with the rank-two tensor q^'^. 
Lowering an index in Ea. (l20p with and using the 
relation 11 “^ = q~^g introduced above, it follows that 
dini) can be rewritten as 

( 22 ) 


This equation puts forward that, like in the f{R, Q) case, 
the metric g^^i, and the metric q^j^i, associated to the con¬ 
nection are related by a matrix H which is fully deter¬ 
mined by the matter sources. It is worth noting that 
defining T = the relation between 5 ^ 1 / and 

q^i, take the form 

q^.. = |t|^(T-l)/5„. ; q>^'' = , (23) 

|T|^ 

which are formally identical to those given in Eq. (fT^ for 
the f{R, Q) theories with T playing the role of S. 

The field equations for the metric can also be writ¬ 
ten in compact form. Starting with the definition of 
written in the form ei?^y(F) = q^i, — 5 ^ 1 ,, we can 
raise an index with the metric q°‘'^ to obtain eR^'^{q) = 
5^" — (fi“^)^ where R^"{q) = i?^a(r)g“‘" is the Ricci 
tensor of the metric q^''. Using now Eq. ([2^ . and the 
fact that = |T| 3 = 2 , we finally get 

^/(9) = [LbiS^’' + T/], (24) 

|T|d-2 


where we have used the fact that the Lagrangian in the 
action (HID can be written as 


Lbi — 


ITI^ - A 


(25) 


Remarkably, the form of the field equations (l24l) for the 
BI theory, which constitutes a well motivated extension 
of GR, is identical with the form found for f{R,Q) the¬ 
ories in m, with the gravity Lagrangian in the f{R, Q) 
case given by L/ = J/2kS. This formal representation 
of the field equations is quite general in Palatini theo¬ 
ries, though some slightly different examples are known 
[sills^. Let us point out that, in the four dimensional 
case it has been observed that BI gravity coupled to 
an electromagnetic field coincides with the action of the 
quadratic theory, not only at the perturbative level (see 
Eq. (|T^ L but to all orders, because the higher-order 
terms cancel out exactly due to the structure of the mat¬ 
ter energy-momentum tensor [^ . 

Like in the /(i?, Q) case. Eg. (Oil) constitutes a system 
of second-order field equations for the metric Given 
that it is algebraically related via Eg. ( 051 ) with g^^, it 
follows that g^y also satisfies second-order equations. It 
also confirms, that this framework naturally accommo¬ 
dates a cosmological constant term, as long as A I. 
In vacuum, = 0, a glance at (( 551 ) confirms that the 
theory boils down to GR plus a cosmological constant 
term, in agreement with this general property of Palatini 
f{R,Q) gravities. 


IV. ELECTROVACUUM SOLUTIONS 

In this section we shall work out black hole solutions 
of Born-Infeld gravity coupled to a spherically symmet- 











ric electric field in arbitrary dimensions, for which exact 
analytical solutions may be found. 
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A. Action and field eqnations 

The matter (Maxwell) action in d dimensions can be 
written as 


5 = J (26) 

where — d^A^ is the field strength tensor of 

the vector potential A^. Let us assume a static space- 
time with a maximally symmetric (d — 2)-dimensional 
subspace (which could be spherical, flat, or hyperbolic) 
with a line element of the form 

ds^ = gttdt^ + gxxdx^ + r(a;)^dn(^_2) , (27) 

where denotes the infinitesimal line element of 

the maximally symmetric subspace and r{x)‘^,gtt, and 
gxx are functions of the spatial coordinate x. In the 
spherical case, dfl(^_ 2 ) = ddf + . 

For static, spherically symmetric, electrovacuum solu¬ 
tions there is a single non-vanishing component, F*^^{x), 
of the field strength tensor. From the Maxwell field equa¬ 
tions, V^F^'^ = 0, this component in the line element 
(HZl) satisfies 


O- = (A + A)^ ; 12+ = (32) 

{X + X)^ 


The functions fl± determine the relative deformation be¬ 
tween and 9 ^ 1 / and play a major role in the properties 
of the corresponding solutions. We now replace Eqs. (ED) 
into the Lagrangian density ([33)l . which becomes 


Lbi = 


d-2 

n+nj -A 




(33) 


From ED, ED and ( 1 ^ the held equations follow im¬ 
mediately after a bit of algebra: 


ei?/(<z) 



These held equations hold irrespective of the number of 
space-time dimensions and arise as an elegant extension 
of the four-dimensional counterparts studied in [^ . 


B. Formal solution of the field equations 

To solve the held equations ([34|) we introduce a line 
element for the metric q^i, of the form 


^ 

r{xY-'^^J-gttgxx 


(28) dY = -A{x)e^^^^^dY + -^dx^ + r^Y)d^l, (35) 


The energy-momentum tensor for an electromagnetic 
held, T/ = - \S/FjFp'^), dehned by 

(|28)l takes the form (from now on n = d — 2) 

r/ = ^ f ■/""" ) , (29) 

OTT y U2xn ^nxn / 

where X = ■ This T/ does not 

depend explicitly on the metric components gtt and gxxi 
which simplihes the resolution of the equations. 

To fully specify the held equations in this scenario we 
need the explicit expression of fl following from Eq. (|2^ , 
which in our case reads 


where '!/f(a;), A(x)^ and f(x) are three functions to be de¬ 
termined by the held equations. Note that f (x) is, in gen¬ 
eral, a function of the radial coordinate x. With the line 
element (1351) we can compute the components of the Ricci 
tensor in d dimensions. The corresponding expressions 
have been worked out in detail in Appendix Equat¬ 
ing such components to the right-hand-side of the held 
equations (IMl) we hnd that Rt — Rx^ = 0 which, from 
Eas. (|A13|) and (IA14|) . implies that Vxx = YxXx, which 
upon integration gives fx = This allows us to in¬ 
troduce a new function A —A/f^ and redehne the x 
coordinate as f'^dx'^ —)> dx^ so that the line element may 
be written in a standard Schwarzschild form 




(A-|-A)/2x 2 6nx2 

02xn (A A)/nxr 


(30) 


ds^ 


—A{x)dY + 


1 

A(x) 


dx^ + x^dVt^ , 


(36) 


where we have introduced the simplifying notation X = 
2 

Given the symmetry of our problem, we now pro¬ 
pose the following ansatz 

h=( ^+^2X2 0„X2 \ ^ (3^^ 

\ 02xn ^—Inxn J 

Since the determinant of is |fi| = we solve (150)1 

to obtain 


where now the role of x as a radial coordinate is apparent. 
This leaves a single independent equation for A(x), which 
from (|A15p and (134)) reads 

I_ [(d _ 3)(A: _ A) - xA,] = , (37) 

where k = —1,0, 4-1 is the curvature of the maximally 
symmetric subspace. In this work we shall focus on the 
case fc = 4-1 for simplicity (see however [3^ for some 
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examples of topological black holes). We now take a mass Putting all these elements into (p 8 )) one arrives at the 
ansatz in d dimensions: following formal expression 


A = l- 


2M(x) 
(d — 


(38) 


where M (x) is the mass function. Replacing this expres¬ 
sion into equation (EZD we find 


A(z) = 1 - 


f l + diG(z) 

I d-3 


(46) 


where 62 is another dimensionless constant given by 


2eM, d- 2 {^--A 
(d-3) \ n_ )' 


(39) 


On the other hand, from the relation ft = g it follows 
that 


dab — 9abf+ : Qmn — ^ ( 4 ^) 

where (a, h) contains the 2 x 2 block and (m, n) the max¬ 
imally symmetric sector (see Appendix for more de¬ 
tails on the notation). It thus follows that the relation 
between the radial coordinate x and the function r'^{x) 
of the metric 5 ^ 1 / is given by 

x^=r^fl_. (41) 


Taking a derivative in this equation, using the relation 

.d — l ~ 

fit = (A -I- A) and gathering terms, we obtain 

^ 

dr \X + X J 

Combining these results with Eg. (15^ we get 


(42) 


2eMr 


„d-2 


- 1 


ft 


1/2 


(A-A). 


(43) 


It is convenient at this point to introduce some simplify¬ 
ing notation. First, since we are mainly interested in the 
case e < 0 (to meet the four-dimensional results stud¬ 
ied in [ 2 ^) we redefine e = where 4 is some length 
scale. The radial variable is rescaled as z = r/vc, where 
^ 2 (ii- 2 ) ^ ^2^2(d-3)^ with = K'^q’^j^Air) represent¬ 

ing a length scale associated to the electric charge. With 
this notation, we get A = — l/z^G-2)^ which is clearly 
dimensionless. With these elements we can integrate (H51) 
as M{z) = Mo{l + SiG{z)), where Mq is a constant iden¬ 
tified with the Schwarzschild mass in vacuum, is a 
dimensionless constant defined as 


_ (d-3)rf-i 

' 2MoZ2 


(44) 


and the function G(z), which accounts for the electric 
field contribution, satisfies 


G. 


M, 

SiMq 


—z 


d -2 


fl_ - 1 



1 


^ 2 {d- 2 ) 


■ (45) 


^ id-S)rt^ 
2Mo 


(47) 


To obtain the physical line element associated to 5 ^^, we 
use Ea.(l40l). which yields 


ds^ 




z^{x)dfl^. 


(48) 


The set of equations (H51) . (HSl and (H51) . with the defini¬ 
tions (liil) and (SZl) and the area function r‘^{x) defined 
implicitly in (1411) . provide a full solution for the problem 
of Born-Infeld gravity coupled to a spherically symmetric 
electric field in arbitrary space-time dimensions. 


C. Exact solution and series expansions 


In order to better understand some important aspects 
of these solutions, before considering the explicit resolu¬ 
tion of (14511 we show how to obtain an explicit expression 
for {x) in dnj. For simplicity we consider an asymptot¬ 
ically flat configuration, without cosmological constant, 
such that A = 1. The 1 case has been considered in 
four dimensions [s^. Using the notation introduced in 
the above discussion, we can write m as 


„2 
' c 



1 

z2G-i) ) 


(49) 


It is clear from this relation that z > 1, i.e., it can not 
become smaller than 1. The equality z = 1 occurs when 
X = 0. This has consequences for Ea. (l4^ . which has a 
divergence at z = 1. This is just a manifestation of the 
fact, well known in four dimensional models, that the 
area function z(x) = r(x)/rc is not a good coordinate 
at X = 0. In fact, from (H^ one can readily verify that 
dr/dx vanishes at x = 0 (or r = Xc), which signals a 
minimum in the function r(x). As is well known, the area 
function r(x) can only be used as a coordinate in those 
intervals in which it is a monotonic function [s^. To 
cover that point, therefore, one should use the coordinate 
X, for which the function G(x) is well behaved. Now, with 
simple manipulations, (I49|) can be written as 


R 

Tc 


{d- 2 ) 


1 


^^2{d-2) 



(50) 


which is just a standard quadratic algebraic equation for 
zd- 2 ^ The solution is straightforward and leads to (see 
















r(x) 

3.0^ 


where 6 d, ad and bd are constants given by 



Figure 1. Representation of r{x) for d = 4 (solid curve), 
d = 6 (dashed curve), and d = 10 (dotted curve). The axes 
are measured in units of Tc. The wormhole throat is located 
at a; = 0, where the area of the (d — 2)-spheres reaches a 
minimum before bouncing into another space-time region. 


Figni) 


“ 2 


(51) 


The modulus in the coordinate x follows from the fact 
that it appears as in the original relation dSD and 
because the range of definition of x can be extended to 
the whole real axis. This gives consistency to our spher¬ 
ically symmetric, sourceless electric field, whose effective 
charge now appears as a topological property of the elec¬ 
tric flux through a wormhole (see also Sec II V El below). 
The throat of this wormhole is located at a; = 0, where 
the area function r{x) reaches its minimum value. 

Now that the domains of the radial coordinate x and 
the area function r{x) have been clarified, we proceed to 
solve dUl) to obtain G{z). We will focus on the asymp¬ 
totically flat case A = 1. Using the binomial expansion 
of the different functions that appear on the right-hand 
side of this equation, an exact solution in terms of infinite 
power series can be obtained for arbitrary dimension. In 
general, the solution takes the form (d > 4) 


Gd{z) = - 


z(3-d) 
zG-i) 


2 F 1 


1 d-3 3d-7 


I 


d- 1 

zG-i) 

d- 1 


2 F 1 


( 1 - 


_d-2’2d-4’2d-4’z2d-4 

I (d-I) d-3 1 ■ 


d-2’ 2d-4’2d-4’ 


^4-2d^(d-l) 


(52) 


where 2^1 is a hypergeometric function. One can verify 
that in the limit z —>■ 00 the GR result is recovered, 
G{z) Ri — (^d'-^_1d+6) other limit of interest is 

when z —^ 1. In that case, we find 


lim Gd(z) 


-^ + ad{z-l)^+bd{z-l)^ (53) 
Od 




ad 


bd 



TT CSC 

TT 

d-2 

d-3 , , 1 

2^{d-2)~^ 

F 

d 

1 

-2 

F 

d—b 

d-2 



(55) 


(2(d-2))^ 
d- 1 


(56) 


where r[a] is Euler’s gamma function. The expansion of 
the metric function (1461) near the wormhole throat be¬ 
comes 


A(z) 


(di-drf) (2(d-2))-^ / {d-3) 

S2Sd (z-I)S V (2(d-2) 


(5i(2(d- 2))”^ 



2 bd{d- 2 ) 

(d-I) 


(z- I)-i"2 


2ad{d - 3) 
{2d-5) 


iz-l) + 


2bd 

(d-I) 


(z- 1 )^^ 


(57) 


One can readily verify that this function diverges at z = I 
if di ^ Sd but is completely smooth if di = Sd- For 
the regular case, the corrections that follow the zeroth- 
order constant term depend on the dimension, being lin¬ 
ear in (z — 1) in d = 4 but going like (z — 1)^^ in 
higher dimensions. Similarly as in d = 4, depending 
on the ratio Si /62 one might find configurations with¬ 
out event horizons which, roughly, can be interpreted as 
two (d—dimensional) Minkowski space-times connected 
through a wormhole. 

We can now consider the form of the metric component 
gtt = —■S^{z)/^+ around the region z = 1, which takes 
the form 


gtt 


(2(d-2))^ 

2 


{Si-6d) (2(d-2))-^ 

S 2 Sd {z — 1)T^ 



aA(2(d-2 ))-^\ 

d2 j 


(Z- I)-i-2 


{5i-5d) (4d-7)^^ ^^^ 


(2(d-2)) 
Sibd (2(d-2))^ 


(d-1) 


(z-I) 


(58) 


For 5i ^ 5d^ the first term in this expression provides the 
main contribution. In terms of the coordinate x, near 
X —>■ 0 it turns into 


(di - Sd) Tc 

2 S 2 Sd \x\ 


(59) 
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For = Sd, we find instead 


1 / addi(2d-4)-^\ /|a;| 

^2 ){r^ 


_ bdSd _ 

d2(d- l)(2(d-2))^ 



(60) 


From this last expansion, it is easy to see that only in 
d = 4 can the near-wormhole geometry be free of cur¬ 
vature divergences, as any derivative of the metric near 
a; —>■ 0 will introduce inverse powers of |a;| even in the 
regular case. Unlike in d = 4, no horizonless solutions 

can exist for d > 4 even if < 1. In that 

case, however, the event horizon is not hiding any interior 
region, as it only exists at the wormhole throat, acting 
in this way as a kind of membrane that separates two 
asymptotically Minkowskian space-times. 


D. Curvature scalars 


To better understand the properties of curvature in 
these geometries, we now compute the Kretschmann 
scalar, K = (see Appendix |B]) around the 

region x = 0 where the deviations with respect to GR 
occur. For arbitrary values of the parameters 61,62 we 
can use Eq. (15^ together with (IBllI) and (IB12I) to expand 
in series around x = 0 {z = 1). We first check the far 
limit behavior, z ^ 1 , which yields Kd ^ in 

agreement with the behaviour of the Kretschmann scalar 
in the RN case of GR. On the other hand, for cc ~ 0, the 
corresponding result admits a generic expression 


^d{6i — 6d)^ Bd{ 6 i, 62 , 6 d) 

B-d — C.9C-9 A,, in + 




r.3d-7 


(61) 


where Ad, Bd, ■ ■ ■ are constants multiplying each term in 
the series. The expansion around x ~ 0 can be easily 
transformed into an expansion in z just by using (EU, 
which behaves there as x ~ (z — The expan¬ 

sion in (EB) shows that the degree of divergence grows 
with d. In the leading-order term we always find a factor 
(di — 5d), which means that the degree of divergence is 
softened in the case di = 6 d- Replacing this choice in the 
metric expressions, and performing again series expan¬ 
sions around x ~ 0 yields 


KdW=5a 


Cdi 62 , 6 d) Ed{ 62 , 6 d) 


^ 2 (d- 2 ) 


^2{d-3) 


(62) 


where Cd, Ed ,... is a set of new constants. This expan¬ 
sion is still divergent for d > 4 (in terms of z one finds 
the same degree of divergence, namely, ~ Cdl{z — 1 )^), 
though it can be further smoothed out if the ratio 62 / 6 d 
is properly chosen, but divergences still persist. How¬ 
ever, for d = 4, the di = 6 d case leads to finiteness of the 


Kretschmann (and other curvature scalars), since one ob¬ 
tains C 4 = E 4 = 0 and 


Ki\si=Sd 




0 {x^) , 


which is a result already found in the literature 


(63) 


M- 


E. Electric flux and action integral 

As pointed out above, from the analysis of the radial 
function r(x) given by Ea. dCT) (see FiglB we concluded 
that our geometry represents a wormhole space-time. We 
can reinforce this interpretation by exploring in some de¬ 
tail the origin of the electric charge associated to our 
solutions. The electric charge can be defined via the flux 
integral 


4* = / *F = uJd-2q , (64) 

where *F is the Hodge dual to the electromagnetic field 
tensor, and the integration takes place over a bound¬ 
ary defined by a (d — 2 ) closed hypersphere with 

cod -2 = 27rT-/r[^^] representing the volume of a unit 
(d — 2) sphere. The hypersurface 5'^“^ contains the min¬ 
imal sphere x = 0 in its interior and, for this reason, the 
net flux is non-zero and yields a finite charge q. The hy¬ 
persphere can be chosen to be on the side x > 0 or on the 
side X < 0. If its orientation is chosen in the direction of 
growth of the (d— 2 )—spheres, then the charge estimated 
by an observer in x > 0 has different sign from that com¬ 
puted by an observer in x < 0. This type of solutions 
were first described by Wheeler and Misner as the result 
of the combination of free electric fields with non-trivial 
topologies [13, giving rise to the notion of mass without 
mass and charge without charge. 

An important property of these solutions is that the 
area of the wormhole, A = ujD- 2 'fc~'^, grows with the in¬ 
tensity of the electric flux in such a way that the electric 
field at the throat is a universal constant independent 
of the particular charge or mass characterizing the solu¬ 
tion. In fact, dividing the electric flux by the area of the 
wormhole throat, we find 


$ 


Wd_2r’c ^ 


1 . 


(65) 


where, recall, h is the scale characterizing the Born-Infeld 
Lagrangian and is proportional to Newton’s constant 
in d dimensions. This puts forward the existence of topo¬ 
logical properties (associated with the wormhole struc¬ 
ture) which are independent of the particular geometrical 
aspects of the solutions. Note that in the limit —)> 0, 
where GR is recovered and the wormhole shrinks to zero 
area, the above quantity diverges. 

Another interesting aspect to note is the fact that the 
spatial integral of the action describing these solutions 
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is finite. Without knowing the explicit form of the so¬ 
lutions, it is possible to evaluate the integral in m 
to obtain a generic result that signals another universal 
property of these solutions. This integral can be written 
explicitly as 


St 








dzz‘^-^{l - X) 
{l+X)T^ 


K^eX 

8tt 

■ {1 + X)T3^ 


Defining where Gd is the d—dimensional 

Newton’s constant and c the speed of light, restoring fac¬ 
tors of Gd and c in (1381) to make it dimensionally con¬ 
sistent (with Mq a mass), and using the definition of 
di, we can rewrite the factor in front of the integrals as 

“ —2Moc^di. Evaluation of the integral gives 
the highly nontrivial result 


St = 2Afoc' 


2 *^1 


dt 


(67) 


where 5d is the coefficient defined in (l54l) resulting from 
the series expansion of the function G{z) in the limit 
z —>■ 1 (or a: —>■ 0). This result is a rather generic prop¬ 
erty of this type of theories and is not restricted to the 
coupling of gravity with Maxwell’s electrodynamics [s^. 
The wormhole structure has been essential to regularize 
the radial integral, which diverges in GR because in that 
case the integration extends down to 0 in its lower range. 
The factor 2 multiplying the mass is due to the need to 
integrate on both sides of the wormhole. Interestingly, 
the resulting action coincides with that of a relativistic 
point particle of mass 2Mq^ at rest. 


V. GEODESICS 

We have seen above in Sec lIVDI that curvature diver¬ 
gences arise generically at the wormhole throat, which are 
frequently regarded as a sign of the existence of space- 
time singularities. However, singular space-times are 
actuallY characterized by the existence of inextendible 
paths [ij, [13 • In particular, if there exist null or time¬ 
like geodesics which cannot be extended to arbitrarily 
large values of their affine parameters, i.e., if they have a 
beginning or an end, then the space-time can be regarded 
as singular. In this section we address this important 
point of the geometries presented above. A comprehen¬ 
sive discussion of the case d = 4 appears in [I^ . 

Taking advantage of the high degree of symmetry of 
our spherical space-times, instead of solving the geodesic 
equation, we work directly with the norm of the tan¬ 
gent vector of geodesics and its conserved quantities. De¬ 
noting the geodesic paths as = (t(r), x(t), ^^(t)), 

conservation of angular momentum implies that we have 
{d — 2) relations of the form Li = gu^- Because of 
spherical symmetry, we can rotate the coordinates so 


that the first (d — 3) angles are fixed at 7 r /2 and all the 
motion takes place in the last angle, 0^-2 = ‘P, so the 
conserved quantity that matters is the total angular mo¬ 
mentum L = r^-^- Time translation invariance gives 
another conserved quantity, namely, E = The 

radial component of the geodesic tangent vector is thus 
( 66 ^overned by 



1 da; 
D+ dr 


lE^ +gtt[ri + 


L 2 

r^(a;) 


( 68 ) 


where 77 = 0,1 for null or time-like geodesics, respectively. 
The ± sign in (ISSl) denotes outgoing/ingoing geodesics. 
From this last equation, we see that the motion in the 
radial direction is analog to that of a particle of energy 
£ = E'^ in a one dimensional potential 


= + (69) 

For null radial geodesics, 77 = 0 = the geodesic 
equation admits an exact analytical expression in the 
form 


± Et{x) 


W[X, d; z] if a; > 0 

, (70) 

2C'[A, d] — W[X, d;z] if a; < 0 


where 
W[X,d-, z] 





d-3 


X 


2^1 



4-2d’ 


1 


(71) 

1 z4-2<i\ 

4-2d’ A J ’ 


and 2 ^ 1 ( 0 , &, c; a:) is a hypergeometric function. In the 
case A = 1, the constant G[X,d] = IF[A, d;l] takes the 
form 


C[l,d] 


3r(4^)r(77^) 


(72) 


A plot of Et{x) is shown in Figl^J from which one readily 
sees that the affine parameter can be smoothly extended 
across the wormhole throat regardless of the parameters 
characterizing the solution. 

We now consider time-like and null geodesics with 
7 ^ 0 and focus on those which can get very near the 
wormhole throat to explore their extendibility through 
the region with curvature divergences. We will thus use 
the approximations given in (15^ and (1501) for gtt in V).// 
(see Eg. (1501) 1. These paths can be classified in three 
groups, depending on the value of di: 


• If di > dd, the potential near the wormhole is an 
infinite barrier, 14// ^ 1/|2;|, and the geodesics do 
not reach the wormhole. In these configurations, 
the wormhole is in a time-like region (below the in¬ 
ner horizon) and, therefore, these geodesics can be 
extended to arbitrary values of their affine param¬ 
eter in much the same way as in GR. 
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E r 



Figure 2. Representation of the affine parameter Et{x) as 
a function of the radial coordinate x for outgoing radial null 
geodesics in d = 4 (solid curve), d = 5 (dashed curve), and d = 
10 (dotted curve). The smooth curves show that radial null 
geodesics can be naturally extended through the wormhole 
throat, X = 0, where curvature scalars generically diverge. 


• If (5i < Sd, the potential is an infinite well, and the 
geodesic reaches the wormhole. Near the wormhole 
gtt oc [see (l59l) ]. and ~ It is possi¬ 

ble to integrate (l68|) in the neighbourhood of the 


wormhole to obtain ± 



'd)) ~ d- 5 / 2 —■ Siiice X g] — oo,-|-oo[, it is clear 
that the range of the proper time t extends over 
the whole real line. 


• If = dd, gtt is either proportional to [xj^"^ or 
to as the wormhole is approached and 14 // 

can be neglected if d > 4 (see [l^ for a discussion 
of the d = 4 case). We thus find that ±E{t — tq) ~ 

^ 2 ^i- 3 ) ’ 'which also allows the affine parameter to 
smoothly go through the x = 0 region. 


We have thus seen that in all those cases in which 
geodesics can reach the wormhole, the affine parameter 
t(x) can always be extended through it to cover its whole 
domain of definition over the real line. We also note 
that the integration of ^ = L/r'^{x) leads to a well- 
defined and continuous (/j(t) because r(x) —>■ Xc > 0 as 
X = 0, being the integrand finite everywhere. Integra¬ 
tion of the coordinate f(r) is also smooth and gives the 
result {t — tc) oc (r — ro)|r — if di < Sd, or 

{t — tc) oc T — To if di = Sdi where 4 is the value of 
the coordinate t at the point of crossing the wormhole. 
These results confirm that the existence of curvature di¬ 
vergences at X = 0 does not prevent in any way the ex- 
tendibility of geodesic paths, which puts forward that 
these space-times are nonsingular. 


VI. CONCLUSIONS 


In this work we have considered metric-affine (or Pala¬ 
tini) extensions of GR including powers of the curvature 
scalar, R, and the Ricci tensor in an arbitrary number of 
space-time dimensions. This scenario generalizes that of 
five-dimensional f{R) theories studied in [l^, and that 
of four-dimensional f(R) [13, f{R,Q) theories [13, and 
Born-Infeld gravities [Ig]. We have explicitly shown that 
the field equations are always second-order and ghost- 
free, as they naturally recover the GR equations in vac¬ 
uum. This property is a result of the independence of 
the connection with respect to the metric. These the¬ 
ories, therefore, provide an appealing new scenario to 
investigate aspects of the AdS/GFT correspondence. 

We have shown that the connection satisfies a set of 
algebraic, non-dynamical equations that result from the 
variation of the action with respect to it. The solution 
to those equations implies that an auxiliary metric 
can be defined such that the independent connection co¬ 
incides with the Christoffel symbols of The field 
equations for can then be cast in Einstein-like form, 
with the right-hand side of the equations representing a 
kind of modified matter source, whose form depends on 
the particular gravity Lagrangian chosen. As the physical 
metric g^v is algebraically related to q^i/ via the matter 
sources, finding a solution for automatically provides 
a solution for 

We have illustrated our analysis employing the gravita¬ 
tional analogue of Born-Infeld theory of non-linear elec¬ 
trodynamics, dubbed Born-Infeld gravity. It has been 
shown that the field equations for this theory admit a 
similar formal representation as those of /(R, Q) theories, 
which confirms the robustness of this formulation. This 
theory naturally accommodates a non-vanishing cosmo¬ 
logical constant as long as A ^ 1. We have obtained 
analytical solutions for Born-Infeld gravity coupled to a 
spherically symmetric electrovacuum field, generalizing 
in this way the Reissner-Nordstrom solution of GR. Re¬ 
markably, we have found that the central point-like singu¬ 
larity of GR is generically replaced by a wormhole struc¬ 
ture. Performing series expansions around the wormhole 
throat, where the modified gravitational dynamics be¬ 
comes important, we have shown that all solutions with 
d > 4 exhibit curvature divergences at the throat. These 
divergences, however, do not prevent null and time-like 
geodesics from being extendible across the wormhole, 
which has a highly nontrivial implication, namely, that 
the resulting space-times are nonsingular. Interestingly, 
though the wormhole size grows with the intensity of the 
electric flux (number of charges), it does in such a way 
that the electric field at the throat is a universal constant. 
On the other hand, the wormhole structure of these solu¬ 
tions plays a key role in the regularization of the action, 
whose spatial integral is finite. 

The results presented in this work, therefore, put for¬ 
ward that black hole singularities in arbitrary dimen¬ 
sions can be cured in a purely classical geometric see- 
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nario governed by second-order field equations. The use 
of metric-affine extensions of Einstein’s theory has been 
fundamental to get this type of modified (but second- 
order) dynamics. The fact that space-times with cur¬ 
vature divergences can be nonsingular is in sharp con¬ 
trast with the widely spread strategy of bounded scalars 
to obtain nonsingular space-times [see e.g. [40l - l45l| and 
references therein], and demands for a reconsideration 
of the role typically attributed in the literature to cur¬ 
vature divergences to characterize space-time singulari¬ 
ties. In our case, the emergence of wormholes connecting 
two copies of the space-time, which give structure to the 
standard GR point-like singularities, has been crucial to 
guarantee geodesic completeness in all cases. We point 
out that the presence of such wormholes seems to be a 
generic prediction of metric-affine geometries, since they 
arise in four-dimensional f{R) and f{R,Q) theories as 
well [H, m, . Further analysis has revealed that the 
wormhole structure in such space-times implies, besides 
geodesic completeness, preservation of causal contact for 
physical observers going through the wormhole [l^. and 
a well behaved (scalar) wave propagation [2^ across the 
divergent region. Such an analysis should be extended 
to the scenario considered in this work. The robust¬ 
ness of these results in similar scenarios with other mat¬ 
ter sources will be explored elsewhere. Extensions of 
these results to models involving further curvature in¬ 
variants, such as powers of the Riemann tensor with 
couplings different from those appearing in the Gauss- 
Bonnet/Lovelock theories, will require the development 
of new methods to deal with the connection equations. 
We hope to address all these issues in future works. 
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where the spherical sector satisfies = 

{x°')rjmn{9'')■ Here the following conventions ap¬ 
ply: a, &, c, d, e are indices corresponding to the 2x2 
sector, while i,j,l,m,n are indices for the spherical 
sector. D is the total number of space-time dimensions 
and n = D — 2 the spherical sector dimensionality. 
The coefficients of the (Levi-Civita) connection for this 
metric read 


J^bam = 0 ; rL„ = -g^^^rdirr^mn ; 

r*, = 0 ; T], = ^5] (A2) 

From the definition of the Riemann tensor we ob¬ 
tain the Ricci tensor 


Rp. = Rpcu = + (A3) 

and we decompose it into its 2 x 2 and spherical sectors 
as follows 


RiT (A4) 

Ri^2 = Rlri - VmnlrOr + {D - 3)g“^VarV,] (A5) 

For the 2x2 sector we consider a static, spherically sym¬ 
metric metric, which can be cast as 


(2x2) ^ { -Ae^^ 0 \ 

9ab 0 A j ’ 

From (IA4I) we obtain for this case 


(A6) 


p(2)_ 

Rtt — 


^2g2^ 


‘^Ax'lj^x 


A 


-^ + 2(^x 


■i’l) 


= -- 
^XX 2 




(A7) 

(AS) 


and therefore it follows that ■ 

On the other hand, concerning the {D — 2) sector, for 
a maximally symmetric space we have R’p^y = k{5‘^rjpi, — 

<5“?7,s^) (with k = 1, 0, —1), which implies = {D — 

Thus, from (IA5I) we have 


Appendix A: Decomposition of the Ricci tensor 

In this appendix we consider the expressions for the 
Ricci tensor corresponding to a metric of the form 


9 


(D) ^ 

fliy 


(2x2) 

9ab 


0 


0 ^ 

(nxn) I > 
ymn j 


(Al) 


R^'^l = [(44 - 3) (fc - V.rVer) - Dr] (A9) 

and raising an index gives Rm^'^ = ^[{D — 3){k — 
g°‘^VarVcf) — Dr]. To make these expressions explicit we 
note that for the metric (IA6I) we have the components of 
the connection + S^^) and P^^, = — ^ 

we find that 
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VtVtr = - 


— + 2'il)x^ 

(AlO) 

2 

Vo;Va;r = Tx 

-^X 

x + ^r. 
Ax 


(All) 

Ur = A 

"^XX 

Tx + i’xTx 

(A12) 


Putting all these elements together we arrive to the equa¬ 
tions 


where C = g°‘^VaTVbf■ Therefore, the Kretschmann 
scalar reads 


+ 4(D - 2) 

2{D-2){D-i){k-g^'^VarVbrf 


(B8) 


Rt — —+ 3—'0X + 2(V'i -I- i^xx) 


(Z?-2)(. 


X -^x 

r A 


2lpr 


Rx^=-^[^+?.^^x+2{^l.l+^xx) 
r, n_5l^ 

■^m — o 

- Ar(^rxx + rx(^-^ 

which are the ones employed in section (IIVBI) . 


{D-m-AvD 

Ax 


(A13) 


(A14) 


(A15) 


Appendix B: The Kretschmann scalar in D 
dimensions 


Consider now the 2x2 sector written as 



—A{x) 0 \ 

0 B{x) ) ’ 


(B9) 


to meet the expression of the metric gab^ corresponding 
to an electromagnetic field, considered in this work. For 
this metric the relevant objects appearing in (IB8I) read 


(V,V6r)(V“VV) = — 


j4| 

4A2 


Bx 


+ 1 ^xx 2B^^ 


Then, the Kretschmann (IB8p reads explicitly: 


(BIO) 


From the components of the connection (IA2I) . asso¬ 
ciated to a metric of the form the non-vanishing 

components of the Riemann tensor @ read 


,a(2) 


R^lf^ = -V.jrg'^^^b^er = -E 


bed ’ Hbj 


p^(^) _ p' 

^bed — 

R^l = Rfl + - SMg'^^VarV^r 

ji(n 


R^ = HSlnVjn - 
R 


a{D) 

ijb 


(Bl) 


*(£’) ^ y VgVbr ^ _^i(D) 


ajb j ^ '-abj 


The Kretschmann scalar is written as 

J^{D) ^ 

Computing each contribution separately one obtains 

(B3) 


(B4) 


pd-D) . _ pj(n) . TX imn(n) 

jmn^H — jmn^H 


- AkC{D-2){D-2,) (B5) 

+ 2C‘^{D -2){D -2,) 

= 2^2(11 - 2){D - 3) (B6) 

= {D- (B7) 


K(^) = k(2) + 

2{D-2){D-3) 

' A 


4A2 “ 


2\2 





(Bll) 


where 


K(2) 


1 

^x 



A 

4^ 

~A 

B + 

ka) ■ 

_2^ 


(B12) 


For the particular case of the line element (l48l) we can 
write A{x) = C{x)/fl+{x) and B{x) = l/(C'(a:)n+(x)). 
Taking into account the relations 12+,^ = H+^rCa,, 

12—,xx — F 12+,r^xx; A^ — (II7./I2+ d2j./12^ 

and Br = —C'r/(C'2l2+) — /{CE^) one can compute 

the Kretschmann scalar (IBllI) , where 


[(712^^,^, -|- Q^Cxx — 12+((70,12+,a; + CQ-^-^xx)] 

(B13) 

and the full expression for is immediately obtained 

by replacing the above expressions into (jBllF 
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